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Very little attention has been paid to the properties of optical wavefronts and caustic surfaces due 
to gravitational lensing. Yet the wavefront-based point of view is natural and provides insights into 
the nature of the caustic surfaces on a gravitationally lensed lightcone. We derive analytically the 
basic equations governing the wavefronts, lightcones, caustics on wavefronts, and caustic surfaces 
on lightcones in the context of weak-field, thin-screen gravitational lensing. These equations are 
all related to the potential of the lens. In the process, we also show that the standard single-plane 
gravitational lensing map extends to a new mapping, which we call a wavefront lensing map. Unlike 
the standard lensing map, the Jacobian matrix of a wavefront lensing map is not symmetric. Our 
formulas are then applied to caustic "surfing." By surfing a caustic surface, a space-borne telescope 
can be fixed on a gravitationally lensed source to obtain an observation of the source at very high 
magnification over an extended time period, revealing structure about the source that could not 
otherwise be resolved. Using our analytical expressions for caustic sheets, we present a scheme 
for surfing a caustic sheet of a lensed source in rectilinear motion. Detailed illustrations are also 
presented of the possible types of wavefronts and caustic sheets due to nonsingular and singular 
elliptical potentials, and singular isothermal spheres, including an example of caustic surfing for a 
singular elliptical potential lens. 
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I. INTRODUCTION 

Among relativistic concepts of direct application to 
gravitational lensing, the observer's past lightcone is per- 
haps the most fundamental. The lightcone concept uni- 
fies the temporal and spatial properties of lensing events 
in a geometrical manner that makes the multiplicity, 
magnification, and time delay of the images arise nat- 
urally. In this view, gravitational lensing by dark and 
luminous matter causes the observer's past lightcone to 
curve into a singular three-dimensional hypersurface that 
self-intersects and folds sharply. It is the singular "fold- 
ing" of the lightcone that is responsible for most features 
of interest in gravitational lensing, such as image mul- 
tiplicity, magnification, and time delay. Additionally, in 
this view, caustics — so relevant to observational lens- 
ing — literally acquire a new dimension, turning into 
caustic sheets (possibly multiple) contained within the 
lightcone itself, and carrying a sense of time as well as 
spatial location. The lightcone concept is most natu- 
rally built by imagining an optical wavefront emanat- 
ing from an observer and traveling out and towards the 
past or future. At the start, the wavefront is convex and 
almost perfectly spherical. Upon encountering gravita- 
tional lenses, the wavefront becomes non-uniformly re- 
tarded, acquiring dents. No matter how small and in- 
significant these dents in the wavefront are soon after 
leaving the neighborhood of the lens, in time the normal 
propagation enhances them, inevitably developing cross- 
ings and sharp ridges. The surface traced out in space- 



time by a wavefront is a (past or future) lightcone. The 
trace of the wavefront 's sharp ridges and "swallowtail" 
points is called a "comoving caustic surface" if thought 
of as embedded in three-space, and a "caustic surface" 
if viewed in spacetime. In other words, a caustic sur- 
face is built out of caustic sheets with cusp ridges and 
higher-order singularities (e.g., swallowtails, elliptic um- 
bilics, and hyperbolic umbilics). Remarkably, in spite of 
the naturalness of this concept, very little attention has 
been paid to the greater significance of the caustics as 
sheets or surfaces. As a result, to our knowledge, the 
properties of optical wavefronts and caustic surfaces in 
gravitational lensing have not progressed much beyond 
the generic local classification of caustic singularities of 
the general theory as developed by Thorn, Arnold, and 
others. Our goal is to initiate a program that investi- 
gates optical wavefronts and caustics as they relate to 
the matter distribution of the gravitational lenses in the 
spacetime. 

The vast majority of gravitational lenses can be mod- 
eled extremely well using a small-angle approximation of 
weak-field perturbations of a Friedmann-Lemaitre uni- 
verse, and projecting the deflectors into planes (thin- 
screen approximation) — see the monographs by Schnei- 
der, Ehlers, and Falco Q and Petters, Levine, and 
Wambsganss M for a detailed treatment. An enormous 
advantage of weak-field, thin-screen gravitational lensing 
is that we can obtain explicit analytical expressions for 
the associated lightcones, optical wavefronts, and caus- 
tics on both the wavefronts and lightcones. Indeed, one 
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of the reasons there has been limited progress in the 
study of quantitative properties of the optical caustic 
surfaces and wavefronts due to gravitating matter is the 
lack of computable physical models with observational 
relevance. Thin-screen, weak-field gravitational lcnsing 
provides us with such a setting. 

Gravitational lensing from a wavefront perspective was 
used by Refsdal || to relate the Hubble constant to the 
time delay between lensed images (see ||, Kayser and 
Refsdal 0, |j| for more). Arnold's singularity theory was 
used by Petters Q to treat the local qualitative properties 
of caustic surfaces (big caustics) in gravitational lensing. 
Friedrich and Stewart 0, Stewart ||, Hasse, Kriele and 
Per lick (9), Low ]l0|] , and Ehlers and Newman [jll) also 
used Arnold's theory to treat wavefronts and caustics in 
general relativity with varied motivations. In the case 
of a strong gravitational field, Rauch and Blandford Q 
numerically computed the caustic sheets due to the Kerr 
metric. Wavefronts in the context of liquid droplet lenses 
have been studied extensively by Berry, Hannay, Nye, 
Upstill, and others — see Berry and Upstill Nye 
and references therein. Readers may also benefit from 



wavefronts 



the popular article by Nityananda |15 
gravitational lensing. 

Our paper calculates explicitly the equations governing 
the wavefronts and caustic surfaces in gravitational lens- 
ing with the standard approximations, and analytically 
relates the wavefronts and caustic surfaces to the gravita- 
tional potential of the lens. This allows us to give a first 
quantitative treatment of caustic surfing, a futuristic no- 
tion suggested by Blandford ]T^] in his millennium essay. 
The motivation for caustic surfing is to lock a satellite on 
a gravitationally lensed moving source so as to observe 
the source at a high magnification over an extended time 
period. We derive the starting equations for the trajec- 
tory to be followed in order to surf the caustic surface. 
The appeal of caustic surfing lies in its potential as a tool 
to obtain information about a distant moving source that 
could not otherwise be resolved. 

Our article is organized as follows. Section II derives 
an explicit expression for the present optical path length 
in units of time of light signals that reach the observer 
in the thin-screen approximation. In Section III, this ex- 
pression is used to define a new lensing map in which the 
sources lie on an instantaneous wavefront, rather than 
on a lens plane. The caustic surfaces and sheets are re- 
lated to the wavefront lensing map in Section IV. Section 
V develops some fundamentals of the concept of caus- 
tic surfing. Applications to the case of nonsingular and 
singular elliptical potential models, as well of singular 
isothermal spheres, are given in Section VI. The latter 
section also gives full classifications of the singularities of 
the wavefronts and caustic surfaces of the noted lenses, 
and explicitly illustrates the caustic surfing concept. We 
conclude in Section VII with some remarks and an out- 
look of future applications. The research reported in this 
paper extends preliminary results on wavefronts in gravi- 
tational lensing first reported by Frittelli and Petters f[% 
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II. THE PRESENT PROPER LENGTH 
FUNCTION AND LENS EQUATION 

We shall calculate the present optical length of gravi- 
tationally lensed light rays from a source to the observer 
using Format's principle and coordinates in space. This 
will include a formula (lensing map) relating the posi- 
tion of a lensed source to the impact positions at the 
gravitational lens of lensed light rays from the source to 
observer. For the convenience of readers who are new 
to gravitational lensing, Sections II. A and II.B.l give an 
introduction overviewing those aspects of lensing rele- 
vant to our wavefront and caustic surface considerations. 
Readers already familiar with the basics of gravitational 
lensing can skip to Section II. B. 2, where we cast the lens 
equation in comoving form. 



A. Cosmological Model for Lensing 

We use the usual cosmological setting of gravita- 
tional lensing, namely, a weak-field perturbation of a 
Friedmann-Lemaitre spacetime. These models are quite 
robust and fit extremely well with observations. A de- 
tailed treatment of the assumptions and approximations 
used to model most gravitational lensing scenarios is 
given in Chapter 3 of B- 

The spacetime geometry in the vicinity of a weak-field 
gravitational lens system (i.e., source, lens(es), observer, 
and light rays between source and observer) is approxi- 
mated by the following spacetime metric: 



= -(l + 2^c 2 dr 2 + a 2 (r)(l-2^^ 
= « 2 (t) 



(1) 



The quantity r is cosmic time, t is conformal time related 
to t by 



dt = 



a(r) 



(2) 



with a(r) a positive function called the expansion factor 
and we use the notation a(t) = a(r(t)). Also <f> repre- 
sents the time-independent Newtonian potential of the 
density perturbation due to the lens(es), and dS^ the 
metric of space with constant curvature K = —1,0,1. 
The weak-field limit is assumed (i.e., |</>|/c 2 <C 1), so 
we ignore terms of order greater than </>/c 2 . Note that 
the potential < fi obey s the cosmological Poisson equation 
— see Section II B 2j . The weak-field assumption near a 
gravitational lens ensures that the bending angles a of 
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light rays deflected by the lens are small, i.e., a <C 1. In 
addition to the latter assumption, we suppose that the 
lens is "thin." This means that along the line of sight, the 
diameter of the matter distribution of the lens or source 
is very small compared to the distances between lens and 
observer, and between lens and source. We shall treat 
the gravitational lens and the point sources as lying on 
planes — called the lens plane and light source plane — 
approximately orthogonal to an axis defined by the ob- 
server and some central point on the lens - referred to 
as the optical axis. All these assumptions allow us to 
consistently restrict lensing to small observation angles 
as measured from the optical axis. 

Now consider the gravitational lensing metric g,Q L in 
(H). There are coordinates (R,6,ip) in space, called co- 
moving coordinates, such that the Riemannian metric 
dS 2 - becomes 



dS 2 K = i _ + R 2 (d6 2 + sin 2 Odtp 2 ). (3) 

We assume that the comoving coordinates are dimen- 
sionless, the expansion factor a(r) has physical dimen- 
sions of length, and the potential </> has physical dimen- 
sions of [length] /[time] . The cosmic time r has di- 
mension of time, while the conformal time t is dimcn- 
sionless. Away from the neighborhood of the lens, the 
gravitational lensing metric g^ L is approximated by the 
Friedmann-Lemaitre (FL) metric 

g£ L = -c 2 dT 2 + a 2 (r) dS 2 ^ = a 2 (t) [-dt 2 + dS^\ . 

This is employs the assumption that the Newtonian po- 
tential cj> decreases to zero sufficiently fast away from the 
lens with the lens-source and observer-lens angular di- 
ameter distances sufficiently large. In a realistic physical 
sense, the Newtonian potential of an isolated body de- 
cays no faster than r _1 . Some plausible ways to justify 
the use of this metric as intended are discussed in |lq| 
and references therein. Still, it is not at all clear how to 
improve on this assumption rigorously. Indeed, there are 
some unresolved technical mathematical issues with the 
standard approximations in gravitational lensing, but it 
is not our intention to address them here. 

In this work we will restrict our discussion to a flat 
cosmological setting. Henceforth, we shall assume that 
the gravitational lensing metric Sq L has K = 0. The as- 
sociated FL metric then reduces to the Einstein-de Sitter 
metric 

gEd = gFL = a 2 (t)dSflat, (4) 

where 

dsl at = -dt 2 +dSl, dS 2 = dR 2 +R 2 {d6 2 +sin 2 Odip 2 ) 

(5) 



with < i? < oo, <6 < 7T, <ip <2ir. Note that for 
K = the comoving coordinates (i?, 9, <p) define spheri- 
cal coordinates in space. The conformal relationship (^) 
yields that the Minkowski metric ds^^ is very important 
for our study, because gEd and ds\ at have the same null 
geodesies up to an affine parameter. The topology of our 
spacetime with a gravitational lens now takes the form 
I X E, where / is an open interval of R and E an open 
subset of R 3 . In most applications of gravitational lens- 
ing, we have E = R 3 — A, where A is a finite set of points 
corresponding to singularities (e.g., point masses) in the 
lens. These singularities would include points where the 
density perturbation of the lens or the Newtonian po- 
tential tj> diverge. This yields a lens plane L that is R 2 
minus a finite set of points. We call E the comoving 
space and r x E the proper three-space at cosmic time 
r. Away from the lens, the metrics on E and r x E are 
the standard Euclidean metric dS$ and metric a 2 (r) dS 2 , 
respectively. When we study the caustic surfaces along 
our past lightcone, we project them into E. 



B. The Present Proper Length Function and the 
Comoving Lens Equation 

1. Physical Setting for Present Proper Length of Light Rays 

We must discuss the basic physical setting needed to 
determine the present optical length of light rays via 
Fermat's principle (see pp. 65-76 of [0] for a detailed 
treatment). Suppose that a light signal emanates from a 
source at cosmic time ts and arrives at the observer at 
the present cosmic time tq. Assume that the observer is 
at the origin of the spherical coordinate system (R, 9, ip) 
in the comoving space S and the source at comoving po- 
sition ps = (Rs,6s,(ps)- Assume that the light ray un- 
dergoes deflection by a gravitational lens on some inter- 
mediate lens plane L. Since the deflection angle at the 
lens plane is assumed small, we approximate the light 
ray's trajectory by an FL piecewise-smooth null geodesic 
consisting of an FL null geodesic from the source to L 
and one from L to the observer. Fix an optical axis 
through an arbitrary comoving point pl on L (e.g., cen- 
ter of mass of the lens); say, pl = {Rl,6l,<Pl) with 
{R,0l,ipl) defining the optical axis. Take p^ to be the 
origin of L. We now consider the set N of all piecewise- 
smooth null FL geodesies that left the source at cosmic 
time Ts and arrive at the observer at roughly the present 
cosmic time to after undergoing deflection at L. The 
travel time differences between these null rays are as- 
sumed to be significantly smaller than the Hubble time 
(1/Hq) and so the scale factor is assumed to change neg- 
ligibly during these time differences. In particular, the 
scale factor is taken to be approximately cl(tq) when the 
different rays arrive. Fixing ts and r , the set Af^°(r) of 
null rays can be parametrized by their comoving impact 
vectors r on L (relative to the origin of L). The rays in 
Af^°(r) will be denoted by vr- 
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But not all the rays in AfJ° (r) represent actual light 
rays. By Fermat's principle (e.g., pp. 66-67, §), the 
actual physical light rays (within our approximations) are 
given by the rays in Af£° (r) whose present proper lengths 
are stationary with respect to variations of the comoving 
impact vector r. To characterize the light rays, we first 
determine a formula for the present proper lengths of the 
rays in M^°{r). Our spacetime with 



gGL = g£ L = a 00 



- 1 



l-2-£]dS b 



. (6) 

can be viewed as an Einstein-de Sitter universe with a 
medium having the following time-independent refractive 
index (e.g., p. 54, @): 



2 

1 



(7) 



Along a null ray v r {s) — (r(s), v r {s)) in Af?°(r), where 
s is a parameter along the ray and i>r{s) is the spatial 
path of vr in E, equation (0) shows that 



n(v r {s)) = 1 - — <j>(t>r{s)) = 



(8) 



where the dimensionless incremental length d£ along 
vr{s) is relative to the Euclidean metric dS$. The right- 
most equality in equation (^) uses the assumption that 
\4>\/ c2 1 an d terms of order higher than </>/c 2 are ig- 
nored. 

The present proper length L TQ (vr) of a ray vr in 
J\f^°(r) is then the ray's optical length relative to the 
refractive medium in the hypersurface tq x E. In other 
words, 



L To (vr) = 



n(i> r ) d£ To . 



Here the increment of length d£ To is with respect to the 
spatial metric a 2 dS$ on tq x E, where 



a = a(t(r )). 

Note that d£ To = ao di. Using ([?]), the present proper 
length of vr in units of time is 



T(v r ) = 



1 -^<t>{pr) ) <U. (<n 



Henceforth, we shall refer to T simply as a present proper 
length function and assume that its values are in units of 
time. 

By Fermat's principle, the light rays are characterized 
by those rays in vr for which T(vr) is a stationary value 
of T. In the next subsection, we shall obtain an explicit 
formula that determines these stationary values. 



2. Present Proper Length Function and Lens Equation via 
Comoving Coordinates 

We shall now provide a practical expression of the 
present proper length (in units of time) in terms of rect- 
angular coordinates in the comoving space E of our lens- 
ing spacetime (/ x E,gcL)- First, let's transform the 
dimensionless spherical coordinates (R, 9, (p) to dimen- 
sionless rectangular coordinates (x,y,z) such that the 
observer is at (0,0,0) and the optical axis (i2, #l,<^j,) 
corresponds to the z-axis. The lens plane L is the xy- 
plane located at z; along the positive z-axis. We shall re- 
label the rectangular coordinates on L as (r, zi), where 
r = (ri,ra) and z\ > is fixed. The source is assumed 
to be located at (s, z), where s — (si, s%). When z > zj, 
the source is behind the lens (and so is gravitationally 
lensed) . 

The cosmological Poisson equation at the present cos- 
mic time r = tq is 



V x 2 </>(*) = 4irGa 2 Q p(x), 



(10) 



where x = (x,y,z), V 2 is the Laplace operator rela- 
tive to x, G is the gravitational constant, and p is the 
cosmological volume mass density. Here p has units of 
[mass] /[length] 3 . Set 

P = a-lp, 4> = a 4>. 
Since (|lCj) is solved by 

cbd 



(x) = -a 2 G/ p(x') 



we have 



(x) = -G p(x') 



|x-x'| 



dx' 



R3 l*-*T 

which has units of [length] 3 /[time] 2 . Integrating jl0| ) 
along the optical axis from the source to the observer, 
we obtain the cosmological 2D-Poisson equation: 



V 2 *(r) = — ^ a(r), 



(11) 



where 



2 f° r° 
^{r) = — / <p(r,z')a Q dz', a(r) = / p(r, z')a Q dz' . 

c J z J z 

Note that \&(r) has units of length, while er(r) has units 
of [mass] /[length] 2 . Let 



*(r) = — / (f)(r,z')dz', a(r) = \ p(r,z')dz' 



It follows that W = \& and a — a^a. Equation ([TJ) then 
yields the comoving 2D-Poisson equation: 



„ 9 ~ , , 8irG _ . , 
V r*M = -5- < r )' 



(12) 
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where and a have units of length and mass, respec- 
tively. The function VP can be expressed formally as a 
solution of ( O ) by 



AG r 

$( r ) = / dr'cr(r') In 

c JR 2 



r — r 



where £o is arL arbitrary, dimensionless, fixed constant. 

Equation (||) yields that the present proper length of 
a null ray v r in A/^S ( r ) separates into two terms: 



o 



T{v r ) = — L(y r ) - #(r) 



(13) 



where 
L(u r ) 



dt. 



*(r) = - 



with L{vr) the (dimensionless) Euclidean length of 
Since bending angles are assumed small, we approximate 
the line integral of <p over i>r by an integral along the 
optical axis, i.e., 



{r,z')dz' 



#(r) 
a 



Note that the scaled comoving potential ^(r) is dimen- 
sionless. 

The contribution of the term aoL(i>r) / c to equation 
( |l3| ) is the time due to the present proper length of the 
light's geometric path relative to the Euclidean metric, 
while the term atfS>(r)/c contributes the relativistic time 
dilation (Shapiro time delay) due to the gravitational po- 
tential d> of the lens. We have 



L(v r )=t x 



(14) 



where ti is the Euclidean length of the segment of Or 
consisting of the straight line from comoving impact vec- 
tor r to the observer and ti s the Euclidean length of the 
straight-line segment from the source to r. Expressing 
the null ray Vr parametrically as vr{s) = (t(s), Or (s)), 
we obtain 



= gEd 



dvr dvr 

ds ' ds 



a\r) 



dOr 



ds 



where the squared magnitude of the comoving velocity 
dvr I ds is relative to d5p . This implies 



ti = 



/•so 


di>r 


ds 






> s L 


ds 


/■to 






/ dt = 


to 









dr 



a(r) ds 



ds 



where tl is the cosmic time when vr arrives at the lens 
plane, so and sl are the respective parameter values cor- 
responding to To and tl, and to = t(ro), ti = t(rt). 
Similarly, 



hs — ti — ts, 

where tg = t(rg) with ts the cosmic time when vr left 
the source. Hence, the present proper length function 
T(yr) can be expressed in terms of conformal time as 



T(v r ) = 2° (t - t s ) - *(r] 
c 



(15) 



We now express the present proper length function in 
terms of our comoving rectangular coordinates. In the 
approximation of small angles, we have 



tl = \/zf + \r\ 2 = zi 



\ r \ 
2^7 



lu = ^/{z-z l ) 2 + \r-s\ 2 = z-z l + 



2[z-ziY 



(16) 



By (|lj) and (|l6|), the present proper length of the null 
ray vr becomes 



T{r,s,z) = T{is r ) = -%T(r,s,z) 



where 



T(r,s,z) = z + 



\jr 

2zi 



2{z - zi) 



*(r). 



(17) 



(18) 



Equation (15) implies that the function T(r,s,z) is the 
conformal time of the null ray vr- Fixing s and z, Fer- 
mat's principle yields that the light rays from source to 
observer are determined by those null rays in J\fJ° (r) with 
comoving impact vectors r satisfying 



grad r T(r, s, z) = 0, 



(19) 



where grad r is the gradient operator in rectangular coor- 
dinates r. By ( |l7| ) and (fi8|), equation (19) is equivalent 
to 



s = — r — (z — Zi)a(r) = s(r, z), 

zi 



(20) 



where a(r) = grad r fy(r) is the change in direction 
through which the light signal vr bends at the lens plane, 
and its magnitude is referred to as the (comoving) bend- 
ing angle. Equation (^p| ) is called the (comoving) lens 
equation. It determines the comoving impact vectors r 
of the light rays from s to the observer. 
For fixed z, equation (^0|) defines a map 

s z : L-> S : s z (r) = s(r, z), 
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called a (comoving) standard lensing map from the lens 
plane L to the light source plane S (i.e., the xy-plane at 
position z on the optical axis), and carries no sense of 
time. This is the comoving version of the cosmological 
lensing map commonly used in gravitational lensing - 
see (28). Note that the Jacobian matrix of s z is symmet- 
ric: 



(Jac s z ){r) = -I- (z-z ; )[Hess r *](r). 



(21) 



Here I is the 2x2 identity matrix and Hessr> & is the 
Hessian matrix of of *S? relative to r = {r\,r2). 

The present proper length of each light ray in Af^°(r) 
depends on the comoving position s of the source on the 
source plane. Inserting j2C| ) in (|l8|), we obtain an expres- 
sion for the present proper length: 



T(r,z) = ^T(r,z), 
c 



(22) 



where 



T(r,*) = T(r,s(r,z),z) 

1 r UI2 



z + 



r / 
z— T + (z - zi)[ \a(r) 
zf V 



r ■ a(r) 



*(r) (23) 



with s(r, z) given by (|2(]). The function T(r, z) has the 
convenient property that the source coordinate z appears 
linearly. We will take advantage of this fact later in this 
work. 

Notice that though equations (|lg|), (|2C|), and (j||) are 
stated for z > zi , they actually hold for all z if we assume 
that i&(r) = when — oo < z < zi. For example, if 
z = 0, then T(r, 0) =0. If a light ray starts out from 
a source at position (s,z) in E with zi < z < oo, then 
r is the light ray's comoving impact vector on the lens 
plane L at position z\. For < z < zi, the vector r 
is the comoving position on L of the light ray's spatial 
path when extended backward from (s, z) to L. If — oo < 
z < 0, then r is the ray's impact position on L when the 
light ray is extended forward from (s, z) to the observer 
to L. Since the case that is of interest to us is when the 
wavefront has passed through the lens, we shall assume, 
unless stated to the contrary, that z > Z\. 



Present Proper Length Function via Angular Diameter 
Distances, Redshifts, and Proper Vectors 



the FL angular diameter distances, redshifts, and proper 
impact vectors commonly used in gravitational lensing. 
Equations @ and (0) yield 



o, 



T(r,8,z) = -Z ili+i, 



*(r) 



a Q e s 



(24) 



Since clq = a(t(ro)) = a(ro) = (1 + zs)a(Ts) and the 
FL angular diameter distance from observer to source is 
ds = a(rs)£ s for K = 0, we get 



a Q t 



1 



d s - 



(25) 



The first term on the right of (EJ) is the time-delay of ir- 
relative to the FL light ray from the source to observer 
in the absence of the lens. Using the formula for the time 
delay (e.g., p. 74 of M and p. 146, H), we obtain 



T(r,s,z) = 



l + z L d L ds 



dh,s 



r 

~d L 

dh^s 
dhds 



s 

*(r) 



1 + zs 



d s 
(26) 



where is the redshift of the lens; c?l is the FL angular 
diameter distance from observer to lens and dL,s the an- 
gular diameter distance from lens to source; r = a(rx)r 
is the proper impact vector of the ray at the lens plane 
at approximately cosmic time T£ and s = airs)s is 
the proper vector of the source relative to the opti- 
cal axis at cosmic time t$\ and ^i,(r) = a^ir) with 
aL = a(t(rr,)). Note that has units of length and 
is the two-dimensional cosmological potential at cosmic 
time r = tl- Fixing s and ds, Fermat's principle yields 
that light rays are determined by grad r T(r, s, ds) = 0. 
Equation ( |26| ) now shows that the latter is equivalent to 
the usual (cosmological) lens equation: 



s = — r - d L S grad r *(r) 
dh 



(27) 



This induces the standard cosmological lensing map in 
the gravitational lensing literature (e.g., [Q, p. 77): 



r, : L - S : r/(r) = -f r - rf L , s grad r *(r), (28) 

where L = {a(T^)r : r 6 L} and S = {a(rs)s : r G S} 
are the cosmological lens plane and light source planes, 
respectively. 



In the majority of the gravitational lensing literature, 
the comoving coordinates of the previous sections have 
not been exploited. In order to make a connection with 
standard gravitational lensing, we now show that the 
present proper length function T(r,s,z) in (|l7|), which 
uses comoving coordinates, can be expressed in terms of 



III. LIGHTCONES, WAVEFRONTS, AND 
WAVEFRONT LENSING MAPS 

Consider a gravitational lensing situation with an ob- 
server receiving light from sources beyond the lens plane 
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where deflectors are located. Even though each light 
source emits perhaps in all directions, generating their 
own lightcones, the observer receives only those rays that 
belong to his past lightcone. Consequently, we shall 
study the observer's past lightcone, rather than a source's 
future lightcone. We are also interested in constant-time 
sections of the observer's past lightcone, which we re- 
fer to as wavefronts. In this section, we shall determine 
parametric equations for the past lightcone and for the 
wavefronts generating the lightcone. 



A. Lensed Lightcones 

The observer's past lightcone C~ in the spacetime 
(/ x E,gGL) is the subset of / x E consisting of all 
past-pointing light rays originating from the observer at 
(tq, 0,0,0). Equivalently, the lightcone Cr is the set of 
all spacetime events in (I x E,gGi>) from which a light 
ray arrives at the observer at cosmic time To. To obtain 
a formula for the points in Cr , consider a light ray vr 
starting out at event (js, s s , z s ), where z s may be positive 
or negative, and arriving at the observer at (r , 0,0,0). 
Here r is the light ray's comoving impact position on the 
lens plane L at Z(, resulting possibly from extending the 
ray backw ards or forwards to L — see discussion at end 
Then the (comoving) lens equations 

show that 



of Section [II B 2 

( pp| ) and present proper length equation 
Vr can be expressed as follows with z as parameter: 



Ur{z) = (r - T(r, z), s(r, z),z), 



(29) 



where z varies from z s to 0, T(r, 0) = 0, and s(r, z s ) = 
s s . Alternatively, equation ( p9] ) holds for a light ray from 
the observer to (ts, s s , z 8 ) by varying z from to z s . 
Hence, by allowing z and r to vary, we obtain the follow- 
ing expression for the observer's past lightcone: 



£~ = {(t - T(r, z), s{r, z), z) : zeR, r E L}. (30) 

The set C~~ is a three-dimensional hypersurface in the 
spacetime (I x S, Qgl) and typically has singularities 
due to distortions caused by the lens. In Section IV, 
we shall give equations that characterize these singulari- 
ties. Notice that in the conformally equivalent spacetime 
(I x 5,g GL ), where 



gGL 



gGL 



= - 1 + 2-hr )dt 



l-2£]dSrf, 



a 2 (t) \ c 
the observer's past lightcone becomes 

L~ = {(T - T(r, z), s{r, z), z) : z e R, r E L}, (31) 

where To = cTo/ao- The observer's future lightcone C + 
and conformal future lightcone L + are given respectively 



as follows: 

£+ = {(t + T (r, z), s(r, z), z) : z E R, r E L}, 
1+ = {(T + T(r, z), s(r, z), z) : z E R, r E L}. 



B. Lensed Wavefronts and Wavefront Lensing 
Maps 

With fixed z, we saw that equation ( pp| ) defines a lens- 
ing map from the lens plane to the light source plane at 
z. There is an alternative interpretation of the lensing 
map stemming from the lightcone point of view. We can 
consider (2C|) and ( |22] ) jointly with the condition of fixed 
conformal travel time, instead of fixed z. 

We define a comoving (optical) wavefront in E as a sur- 
face of constant conformal present proper length T(r, z). 
One should keep in mind that conformal time flows at a 
rate different from that of cosmological time — see equa- 
tion (^|). However, the surfaces of constant conformal 
time T are the same as those of constant present proper 
length T; they are only labeled differently — see (|22j). 
Since we are not particularly interested in the specific 
labels of the wavefronts, which are not observable, we 
prefer to use the conformal time 

T = T(r,z) 

to describe the wavefronts. We shall refer to T as a con- 
formal present proper length function. 

Due to the linearity of (E3h in the variable z, we obtain 
a parametric expression in terms of (r, T) for all the z- 
coordinates of the comoving wavefront by solving for z: 



T - r-a(r) + — \a(r)\ 2 + *(r) 



z(r,T) 



\a(r)\ 



r ■ a(r) 



(32) 



2zf 



zi 



The values of z(r, T) are assumed to obey z(r,T) > zi- 
In other words, we are interested primarily in the lensed 
portion of the front, i.e., the portion of the front with 
z-values on the optical axis from the lens plane onward. 
Note that if we allow z(r,T) < zi, then it is assumed 
that ee in @ — see the end of Section |HBl . 

Equation ( p2| ) can now be used in (|2^) to obtain an 
expression for the points of the wavefront that are on the 
xy-plane at position z(r,T): 

s(r,T) ee s(r,z(r,T)) 



= z(r,T)y- - a{r)j + z t a{r). (33) 

In other words, the pair (s, z) determined by ( |33"| ) and 
( |32| ) defines a chart (possibly with singularities) on the 
portion of the wavefront that went through the lens 
plane. Of course, the pair can also be used as a chart 
for the portion of the front behind the lens plane (e.g., 
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by setting ^'(t*) = 0). The comoving wavefront at con- 
formal time T is then given by 



W(T) = {(a(r, T), z(r, T)) : r e L} C £, 



(34) 



which is a two-dimensional surface that generically has 
singularities due to lensing. Notice that if we allow — oo < 
z < oo and < T < T , then W(T) traces out the 
conformal light cone in (j3l|). 

For the cosmological present proper length 



t = T(r,z), 

the wavefront is given in the spacetime (/ x E,^gl) at 
cosmic time tq by 



W (t) = {(r - t, s{r, t), z(r, t)) : r e L} C / x E, 

where s(r,t) = s(r,T(t)) and z(r,t) = z(r, T(t)) with 
T(t) = ct/a Q . Equation ( po[ ) shows that Wo(t) is a sub- 
set of the past lightcone C~ . The wavefront Wo(t) is 
the set of positions of all sources whose light rays left at 
the same cosmic time To — t and arrive at the observer 
at roughly the present cosmic time To. Most observed 
lensed light sources are at distances ranging from of or- 
der tens of kiloparsecs (e.g., microlensing) to Gigaparsecs 
(e.g., quasar lensing). The observation time is typically 
of order weeks to a few years. Consequently, we expect 
the wavefront Wo(t) to change negligibly during the ob- 
serving time, unless the front happens to pass through a 
higher order singularity during the observation. 

The wavefront Wo(t) lies in the proper space {tq — 
t} x E. Projecting the front into the comoving space, we 
obtain 



]< it) = {(a(r, t), z(r, t)) : r e L} C E, 

which is a constant cosmic time slice of C~ (namely, To — 
t = constant). The wavefront W§(t) differs from W(T) 
in E merely by the label T = ct/ao- The two wavcfronts 
are isometric as singular spaces in the comoving space 
with the Euclidean metric. 

Now for fixed T, the lens equation (^) determines a 
mapping 

w T : L -> W(T) C R 3 , 
called a comoving wavefront lensing map, defined by 



w T (r) = (s(r,T),z(r,T)). 

By allowing for z- values with z(r,T) < Zi, the map wt 
can be separated into two single-valued maps 



w+{r) i£z(r,T)>z h 
w^{r) if z(r,T) < z u 



where w+ : L -> W+(T) and w T :L^> W~(T) with 
W+(T) the subset of W(T) for which z > z h while 
W _ (T) is the subset with z < Zi. In other words, the map 
Wt defines two (possibly singular) coordinate patches on 
W(T). Since W~(T) did not pass through the lens, it 
has no singularities and so the map w^, is of little inter- 
est for our purposes. Instead, we shall focus primarily on 
wt- Unless stated to the contrary, we shall assume that 
wt{t) — w^p{r). 

It is important to emphasize that rather than the usual 
mapping from a lens plane to the light source plane (such 
as given by the comoving and cosmological lensing maps 
s z : L — » S and rj : L — > S), we now have a new mapping 
- one from the lens plane to a wavefront. This gives 
us an alternative interpretation of gravitational lensing 
where the light source plane is dispensed with and sub- 
stituted with a wavefront: the locus of points that can be 
reached in a given time by light signals emitted simulta- 
neously from the observer in all possible directions, after 
passing through the lens plane. It should also be kept 
in mind that even though the wavefront as a surface is 
different from a source plane, in typical lensing scenar- 
ios the wavefront surface lies very close to a plane in the 
region of interest, namely around the optical axis. 



IV. CAUSTICS ON WAVEFRONTS AND 
LIGHTCONES 

As a surface in the comoving space E C R 3 , the wave- 
front W + (T) typically develops singularities — called 
caustics — due to the distortions caused by the lens. 
More precisely, a point r in the lens plane L is a critical 
point of the wavefront lensing map 



wt ■ L — > W+(T) C R 3 



w T {r) = (s(r,T),z(r,T)) 



if rank[(Jac iotX?")] < 2, where Jaci^T is the Jacobian 
matrix of wt and z(r, T) > z\. The set of critical points 
of wt will be denoted by Crit(wr). The set of caustic 
points of wt or on W + (T) is the set 

Caustics[W+(T)] = w T [Crit(w T )} 

of all critical values of wt- 

The Jacobian matrix of wt is the following 3x2 ma- 
trix: 



Jac wt 




(35) 



where z = z(r,T) is given by ( p2| ) and s = s(r,T) by 
). Define a function / : L — » R by 



fir) =zi + ^ - *W. 
2zi 



(36) 



9 



Note that f(r) — £i — ^f(r) within our approximations 
(see (|l6|)), so ao/(r)/c is the present proper length of 
a light ray from the point r on the lens plane to the 
observer. By (po|), we see that 

s(r,T) = s(r,z(r,T)) = z(r,T)V r f(r) + Zl a(r) 

Consequently, the Jacobian matrix of s(r, T) relative to 
r can then be expressed as 



ds 
dr 



ft fe + zfll + Z&u h WI + Z hl + *i*12 



(37) 

where /j, /y, and $y denote the usual partial derivatives 
relative to r = (ri, T2)- Note that the Jacobian matrix in 
( |37| ) is not symmetric, unlike the usual Jacobian matri- 
ces ds z /dr and dr//dr of the comoving and cosmological 
lensing maps. If z is fixed, then (|37j ) reduces to the usual 
comoving symmetric case — see equation (^l[): 



ds(r,T) ds z {r) 



Or 



dr 



zi 



-I - (z - z{) (Hessr 



We have rank[(Jacu>x)(f*)] < 2 if and only if every 
2-squarc minor of Jacti>T vanishes: 



dct 



ds 
dr 



dot 



dct 



dr 

dz 
dr 
ds 2 
dr 

dz 
dr 



dsi dz 
dr± dr 2 



ds 2 dz 
dr\ dr 2 



dsi dz_ _ 
dr 2 dr\ 



(38) 



(39) 



ds 2 dz 
dr 2 dr\ 



= 0, (40) 



where s — (si,s 2 ). Since the conformal time T is fixed 
when studying W + (T), we have 



T(r, s(r, T), z(r, T)) = T = constant, 
which yields 



9T 
dr 



dT ds 
ds dr 



dT dz 
dz dr 



0. 



By (19), we get 



dr ' 



(41) 



(42) 



while equation ( fl8| ) yields dl/ds = (s — r)/(z — z{). The 
latter vanishes if and only if a(r) = vjzi (apply ([42"|) or 



(po|)), which holds if and only if r is a critical point of 
the function / in (j36|). Note that for r to be a critical 
point of £ it must solve 

r 9/ 

— - a(r) = ^-(r) = 0. 

Generically, the critical points of / are nondegenerate 
and so are isolated points (e.g., p. 240, Petters et al.). 
For this reason, we shall assume — unless stated to the 
contrary — that r ^ s. Consequently, 



dT 
7Ts 



^0, 



z - zi 
Equation ( pl| ) yields 

dz 
dr 



dT 

dz 



2( Z - Zl y 



* 0. (43) 



(gT/gg) 

(»r/a«) 



9s 
<9r ' 



(44) 



where by ( p3[ ) the bracketed term is nonzero. Plugging 
dz/dri and dz/dr 2 from (|44|) into ( ^9|) and ([h]) implies 
that 



ds i 9z 
9ri dr2 
ds2 dz 
dr\ dr 2 



dsi dz 
dr 2 dri 
ds 2 dz 
dr 2 dri 



(dl/ds 2 ) 
(dT/dz) 
(dT/d Sl ) 



det— , 
dr 



{dT/dz) 



dct ■ 



9r 



(45) 



(46) 



Equations (||) - (gcj), @, @, along with our assump- 
tion in (^3|) , then show that the Jacobian matrix of wt 
has rank below 2 if and only if the Jacobian determinant 
of s vanishes. Thus, at the conformal time T the caustics 
on the wavefront are given by 



ics[W(T)] = {(s(r,T),z(r,T))}, (47) 



Caust 



where r G L, s(r,T) is given by ( p3| ) with T fixed, and 
s = s(r, T) is subjected to the constraint 



detg(r,T) = 0. 



(48) 



Equation ([47]) extends to wavefront lensing maps the 
usual concept of caustics in gravitational lensing, where 
lensing maps are from a lens plane to a light source plane. 
Allowing the conformal time T to vary in (^) yields 
evolving caustics on the wavefront as the front propa- 
gates. The caustics on a wavefront trace out a caustic 
surface on the conformal past lightcone L~ is given by 

Caustics[L-] = {(T -T(r),s(r,T(r)),z(r,T(r)))} 

_ (49) 

where To = ctq/oq, r e L, and s(r,T) obeys (Hq). Since 
T varies in the case of a caustic surface, the vanishing- 
determinant condition (Esh forces T to be a function of 
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r, which we have expressed by T(r). Also, note that the 
caustics lie on the portion of with z > z\. Projecting 
Caustics[L _ ] into the comoving space E, we obtain the 
comoving caustic surface 

Caustics*^"] = {(s(r,T(r)),z(r,T(r)))\, (50) 

where rgL and s(r,T(r)) satisfies (f48|). In the space- 
time (/ x E, gGi), the caustic is given by 

Caustics[/T] = |(t - r(r), s(r,T(r)), z(r,T(r)))\ 

where tq is the present cosmic time and r(r) = 
aoT(r)/c, r e L, and s(r,T(r)) obeys (|||). Projecting 
Caustics [C~] into E yields a surface that is isometric to 
CausticssfL - ] as singular spaces relative to the Euclidean 
metric on E. Analogous to equations ( |49"1) — (H^), the 
caustics of the future lightcone and their projection into 
the comoving space are given as follows: 



Caustics LC 



+1 - 



{(7b+T(r),a(r,T(r)),z(r,T(r)))}, 

(52) 



Caustics[L+] = {(T + T(r),8(r,T(r)),z(r,T(r)))}, 

(53) 



Caustics £ [L+] = {(a(r,T(r)),«(r,T(r)))}, 



(54) 



where reL and s(r,T) satisfies (p8|). 

The previous discussion traces out the caustic surfaces 
of an observer's past lightcone using constant time slices. 
The caustic surfaces can also be traced out using constant 
z-slices. Slices of the caustic sheet by constant z-planes 
are curves on the light source plane, and coincide with 
what are commonly referred to in the lensing literature 
as "caustics." We shall refer to such caustics as z-sliced 
caustic curves. These caustic curves can also have sin- 
gularities, such as cusps. It is important to add that our 
discussion shows that points on a z-sliced caustic curve 
actually occur at different cosmic (or conformal) times. 
Often times, points on caustic curves are treated in the 
lensing literature as if they occur at the same cosmic 
time. 

In order to calculate the points on the comoving caustic 
sheet, we impose the condition that the Jacobian of the 
lens map be singular — see ([[o]) and (|54|): 

det^=0. 
or 

This condition can be imposed either at constant z or 
constant T if the intention is to produce the caustic sheet. 
For constant T, we have that s is given by the wavefront 
lensing map (|33|), i.e., s — s(r,z(r,T)), while for con- 
stant z the function s is given by the standard lensing 



map (J20|), i.e., s = s z (r). In the latter case, we consider 
the caustics of the mapping s = s z from the lens plane 
to the light source plane, with the intention of eventu- 
ally letting the light source plane at z sweep (along the 
optical axis) the entire range behind the lens. 

By equation ([H]) — equivalently, equation ( |37j ) with z 
fixed — we see that the Jacobian matrix of s z depends 
linearly on z: 



ds z 
dr 



l + (z- Zl )h 



-(Z - Z;)*21 



l+(z-z t )f. 



22 



(55) 



where / is given by ( p6[) . Consequently, the determinant 
of ( |55|) is a quadratic in z. Explicitly, the vanishing of 
the Jacobian determinant of s z is 

= l + (z- zi^zf 1 - #ii - # 22 ) + (z - ztfUz- 1 



-*h)(z, 1 - * 22 ) - *ij 



(56) 



There are two solutions for z as a function of r, repre- 
senting the z~ coordinate of points on the caustic surface: 



z±(r) = zi + 



*H + * 22 - 2Z^ ± J(* n - * 22 ) 2 + 4f 2 2 



2[(z r 1 -* 11 )(z r 1 -* 22 )-^ 2 ] 



(57) 



Evaluating s z — s(r,z) at z±(r), we obtain the trans- 
verse coordinates of points on the caustic surface: 



s±{r) = s{r,z±(r)) 



z±(r) 



r - \z±(r) - zi) a(r). 

(58) 

By varying r across the lens plane, the pair 
i s ±( r )7 z ±( r )) traces out the comoving caustic surface 
for either the future or past conformal lightcone: 

Caustics£[L + ] = |(s±(r), z±(r)) : r € Lj 

— CausticssfL - ]. (59) 

V. SURFING A CAUSTIC SHEET 

In his Millennium Essay, Blandford speculated 
about some possible novel gravitational lensing ways of 
probing the cosmos. One of these dealt with caustic 
sheets: "Our observations need not be passive. .. For 
example, suppose that we launch an array of robotic tele- 
scopes and use three of them to measure the velocity of a 
caustic sheet (from by a bright source) as it passes Earth; 
a fourth telescope could be made to "surf" the wave and 
observe the source with considerable magnification for a 
long time." 



11 



We now apply the results of the previous section to 
show how a telescope may surf a sheet of a caustic surface. 
This will be done in two steps. 

First, change the point of view by assuming that the 
telescope is the observer. We shall then consider the past 
comoving caustic surface of the telescope and the future 
comoving caustic surface of the source simultaneously, 
and determine the analytical form of the equations for 
these surfaces. Assume that the telescope lies somewhere 
on the future comoving caustic surface of the source (i.e., 
the source is seen at an extremely high magnification 
corresponding to being on a caustic). For the time being, 
we shall also suppose that the source is at rest relative 
to the lens. The past comoving caustic surface of the 
telescope is given by an expression of the form 

5 (x) = 0, 

where x = (s, z). In principle, one can find So by solving 
for r in terms of s± using ( |58| ) and inserting those values 
of r into ( |57| ) to obtain z± in terms of s±, say, z± = 
F±(s±). In this case, the past comoving caustic surface 
of the telescope is given by 



S (x) = z ± - F ± (s ± ) 



0. 



(60) 



Since we are using the small-angle approximation, we can 
approximate the future comoving caustic surface of the 
source by (pQ) if s and z are transformed as follows: 



s — > — s 



where £ti and l s i are the telescope-lens and lens-source 
Euclidean distances, resp. Hence, we suppose that there 
is an expression for the future comoving caustic surface 
of the source of the form 



fib(x) = 0. 



(61) 



Second, suppose that a source moving relative to the 
lens with slowly varying 3-velocity V in comoving coordi- 
nates is observed extremely magnified by a telescope at 
a given instant of conformal time, say, t = 0. Assuming 
that the telescope lies somewhere on a caustic sheet of the 
source away from singularities (e.g., cusp ridges, swallow- 
tails, elliptic umbilics, hyperbolic umbilics), what should 
the subsequent motion of the telescope be in order to 
track the image of the source at a peak brightness? Sup- 
pose that during a time increment St the source moves 
keeping the distance to the lens plane approximately the 
same, and assuming its motion to be relatively slow, the 
caustic sheet of the source in motion will not differ from 
that of the source at rest other than by a general trans- 
lation in the direction of motion, i.e., 

S(x,St) = S (x-V St) = 0. 

If the telescope lies on the caustic sheet at point Xo at 
time t = 0, then 

S(x ,0) = = S o (x o ), 



For the telescope to see a bright image of the source dur- 
ing a length of time St, the telescope needs to move with 
velocity v to another point on the caustic sheet, i.e., to a 
point x = Xo + v<5t so that 

S(x,<5t) = = S (xo + (v-V)5t). 
Taylor expanding to first order in St, this is equivalent to 



VS • (v - V) = 0, 



(62) 



where VSo is evaluated at the original position of the 
telescope Xo and • stands for the Euclidean scalar prod- 
uct of the two vectors. This means that the telescope's 
velocity must differ from that of the source at most by a 
vector tangent to the caustic sheet. Clearly v = V would 
be one way to stay on the caustic sheet, but it may require 
more energy than necessary. We want the smallest speed 
needed to stay on the caustic sheet. In order words, we 
shall minimize |v| 2 for a fixed spatial position and time 
subjected to the constraint in (p2[). This will be valid 
for the time increment St, i.e., we are considering only 
the zeroth iterate of caustic surfing. The vanishing gra- 



dient of the Lagrangian C 



■A VSo-(v-V) 



with respect to v = (v±, i>2,«3) and the multiplier A yields 
A = — V • VSo/|VSo| 2 . Hence, we obtain the minimum 
velocity 



-AVS 



V- VS 

WW 



VS , 



(63) 



where VSo is evaluated at xo. In other words, the min- 
imum velocity is the projection of V to the normal vec- 
tor VSo to the caustic surface at xo- This is physically 
the minimum speed since VSo points along the shortest 
direction between the caustic sheet at time t = and 
the sheet at t = St. The speed |v TO j„| might be consid- 
erably smaller than the source's speed |V| depending on 
the circumstances. The normal unit vector to the caustic 
VSo /I VSo I can be calculated either from the expression 
d6l| ) if available or from the parametric version of the 
caustic map given by ( |57| ) and (58). 

In a physically realistic situation the telescope will 
spend some time to reach this velocity, and one will need 
to adjust the velocity continuously even if the source is 
moving at constant speed, due to the curvature of the 
caustic sheet. 

We shall illustrate the caustic surfing concept with a 
singular-elliptical-potential example in Subsection VI B 



VI. EXAMPLES AND ILLUSTRATIONS 

In this section we illustrate explicitly the construction 
of the wavefronts in some of the most widely used models 
for the deflection potential in astrophysics. In all the 
cases that we deal with the expressions for the deflection 
potentials stated are assumed to hold only in the vicinity 
of the optical axis. This assumption is maintained in the 
construction of all the figures. 
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A. Nonsingular Elliptical Potential Lens 

In this subsection we illustrate the case of an elliptical 
potential, which has the following form: 



i> ep (r) = A y / r2 + (l-e)r2 + (l + e)r 



where Aq is a dimensionless constant. This potential is 
often used to model elliptical galaxies with Aq propor- 
tional to the velocity dispersion of the lens, whereas the 
dimensionless parameters r c > and e > determine the 
core radius and ellipticity of the lens. A non-vanishing 
core radius r c guarantees a non-singular surface mass 
density, which in turn ensures that no light rays are ob- 
structed in their passage through the lens. The bending 
angle is given by 



a ep (r) = Vrip(r) 



A ((l + e)n, (l-e)r 2 ) 
^Jr 2 c + (1 - e)r\ + (1 + e)r\ 



For fixed parameter values of Aq, r c , and e, these two 
explicit expressions for the potential and bending angle 
can be used in z(r, T) and s{r,T) to produce plots of the 
constant-time wavefronts embedded in three-dimensional 
space beyond the lens plane. At early times beyond the 
lens plane, the wavefront is essentially spherical, but as 
it progresses away from the lens plane the central sec- 
tion of the wavefront begins to lag behind the rest of 
the wavefront, eventually folding and developing multiple 
sheets. The progression of the wavefront in time shows 
different regimes, from smooth to various singular types. 
In the following, we produce a representative surface of 
each regime, effectively classifying the singularities of the 
wavefronts of an elliptical lens. This extends the brief 
summary of such wavefronts given in 0. Note that for- 
mally our treatment is along the future lightcone of the 
observer, which can be interpreted as the past lightcone 
in a time coordinate running towards the past, because 
the spacetime is static. 



1. Wavefront Singularities 

We arbitrarily fix the values of the parameters at 
A Q = 1.01, zi = 90, e = 0.0002 and r c = 1.001. These val- 
ues are chosen for purely pedagogical reasons. Notwith- 
standing, in making our choice we took care to ensure 
that the approximation of small angles was met, in order 
for the plots to be qualitatively representative of lensing 
problems. 

Let Ti be the time at which the ray that passes through 
the center of the lens (i.e., the origin) reaches the lens 
plane. With our choice of parameters, it takes the value 
Ti = 90. We find that for times at least as late as T = 
Ti + 8900, the wavefronts are regular on the other side of 
the lens plane, but sometime before T = T; + 9800 the 



first singularity occurs in the form of a single point on the 
wavefront. The time scales in this particular example are 
irrelevant; in physically accurate lensing situations the 
time scale should agree with the distance scale to the lens 
plane in order of magnitude. The time of the first singular 
wavefront is about T = Ti + 9699.27, and is calculated in 



(64) Section VI A 2. After the first singularity, the wavefront 



successively goes through four more distinct singularity 
regimes. 

The regular regime is illustrated in Figure |l|, which 
shows a field view of the wavefront. The observer is lo- 
cated to the far right, at (0,0,0) in the plot coordinates. 
The wavefront is distorted with respect to a sphere, and 
develops a dent that points toward the observer along the 
optical axis. The scale of the optical axis is greatly mag- 
nified (about 10 4 times) with respect to the transversal 
axes in order to better show the retardation effect on the 
wavefront. The "spike" has a total length of about 10 -3 
on a sphere of radius about 9 x 10 3 , so the effect is small 
in the proper scales and is well within the approximation 
of small angles. The tip of the "spike" is perfectly smooth 
and presents a convex surface to the observer. There are 
no singularities anywhere in the wavefront at this time 
(T = 7} + 8900). 

At a later time, the first singular regime develops, il- 
lustrated in Figure E|. The tip of the spike, which ear- 
lier was convex toward the observer, develops a self- 
intersection and a sharp cuspidal ridge in the form of 
horizontal "lips." Our use of the term "lips" in this con- 
text is entirely for descriptive purposes and should not be 
understood in the technical sense used in Arnold's the- 
ory, in which the term "lips" refers to the As(+) planar 
caustic metamorphosis (e.g., pp. 375-376, 381 of B). 
The surface of the wavefront that faces the observer is 
now concave, and limited by the two cuspidal ridges out- 
lining the "lips." The two cuspidal ridges pinch off the 
back convex sheet of the wavefront at two symmetrical 
swallowtail points. The swallowtails are best appreciated 
in the top panel, where a side view of the wavefront is 
shown. The bottom panel shows a front view, where the 
lip character of the cuspidal ridge can be appreciated. 
The wavefronts in this regime all have this horizontal lip 
singularity, but the "lips" grow in time continuously after 
starting out from a single point at the center. The pic- 
ture shows the wavefront at T = TJ + 9800, which is still 
quite early in this regime. One notable fact about this 
wavefront is that, locally, it is indistinguishable from the 
early singular wavefronts in the implosion of a triaxial el- 
lipsoid [|18|, in which case, the cuspidal lip-ridge lies "in 
the outside." For the gravitational lensing case that we 
are illustrating, in a global sense, the cuspidal lip-ridge 
lies "inside." 

The next singular regime is illustrated in Figure |§|. The 
point at the center of the concave sheet of the wavefront 
facing the observer turns momentarily singular at ap- 
proximately T = Ti + 10148.84, and grows another lips 
singularity, this time vertically oriented. These vertical 
"lips" ride on top of the now mixed convex-concave sec- 
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tion of the wavefront and present a concave surface to 
the observer, as can be appreciated from the two slices in 
the middle and bottom panels. The middle panel shows 
a vertical slice of the wavefront through the optical axis. 
It may be difficult to appreciate that the foremost sheet 
pinches off the underlying originally existing sheet well 
before touching the originally existing cuspidal ridges at 
the top and bottom. The bottom panel shows an optical- 
axis horizontal slice, where the cuspidal character of the 
newly created forefront sheet is evident. The two ends of 
the vertical "lips" are two symmetrical swallowtails. All 
the wavefronts in this regime have two perpendicular lips, 
yielding a total of four swallowtails. The perpendicular 
lips grow continuously in time, but the inner lip starts 
at the center and grows until it touches the outer lips. 
This picture shows the wavefront at T = 10200, which is 
quite late within this regime since the inner lip is about 
to merge with the outer one. 

The next regime consists entirely of a single wavefront 
at the critical time when the inner lip singularity merges 
with the outer one. This wavefront is represented in Fig- 
ure ||. The singular ridge on the wavefront has the general 
outline of a football in a vertical position, inscribed in an 
outer astroid or diamond. The "football" outline is the 
distorted last expression of the inner lip-ridge, whereas 
the astroid is the distorted last expression of the outer 
(original) horizontal lips. At this time the two cuspi- 
dal ridges merge right before they undergo a permanent 
change. The wavefront at this critical time has two swal- 
lowtails symmetrically located along the horizontal di- 
rection, and two singular points symmetrically located 
along the vertical direction where the two lips merge. In 
the context of wavefronts, these points have no greater 
significance than the merger of a cusp ridge and a swal- 
lowtail surface. However, their greater significance is that 
they are indicators of the presence of hyperbolic umbilic 
points in the caustic sheet, which we describe in the next 
subsubsection. The figure represents the time T = 10420, 
which is close enough to the critical time at the resolution 
that we are using. 

The last regime comprises all the wavefronts later than 
the critical time at which the inner "lips" merge with 
the outer "lips." All such late wavefronts have a "regu- 
lar" cuspidal ridge delimiting a concave surface that faces 
the observer, namely, a cuspidal ridge which is smooth 
as a curve in three space. Behind the foremost concave 
cap, the wavefront self-intersects and has another cusp- 
idal ridge with four singular points in the general shape 
of an astroid or diamond. One such wavefront at time 
T = T[ + 10900 is shown in Figure |. In this front view, 
the foremost oval cuspidal ridge is clearly distinguish- 
able, as are the two swallowtails in the back. However, 
there are two other swallowtails completing the diamond, 
which lie behind the foremost oval cap and are blocked 
from view in this picture. The qualitative structure of the 
wavefront in the diamond-ridge neighborhood is shown in 
the bottom panel of Figure pi 

A very late wavefront at time T = T; + 55000 is shown 



in Figure |6|. It can be seen that as time moves on, the 
oval and diamond cuspidal ridges remain essentially un- 
changed except that the oval ridge in the foreground 
grows at a much faster rate than the diamond-shaped 
ridge in the back, and the wavefront eventually acquires 
a shape resembling that of a goblet. The figure shows 
both a field view of the "goblet" and a magnification of 
the goblet's throat to show the diamond ridge. The di- 
amond ridge structure is almost ubiquitous in wavefront 
evolutions that are not axi-symmetric. For another view 
of one of these local diamond ridges in wavefronts, we 
refer the reader to the bottom panel of Figure [h]. The 
global "goblet" wavefront is typical of spherically sym- 
metric regular potentials, with the major difference that 
the throat of the goblets in the spherically symmetric case 
degenerates down to a single point, lacking the compli- 
cated diamond structure for the case of elliptical sym- 
metry. The diamond ridge represents the nondegenerate 
version of the throat of the goblet. The goblet's throat 
has three wavefront sheets, each generically giving rise to 
an image of the light source. 

What we have described in this section are precisely 
the singularities and metamorphosis of wavefronts in the 
observer's lightcone in the case of an elliptical poten- 
tial. In this case, we found that the wavefronts have 
three types of singularities: cuspidal ridges, swallowtail 
points, and points of transversal self-intersections. These 
three singularity types are shared by generic wavefronts 
in space — see p. 55 of [|l9|. We also showed that two 
kinds of metamorphosis occur for the wavefronts due to 
an elliptical potential. In the first place, we have found 
two occurrences of the birth of "lips" on the wavefront. 
Secondly, we have found two hyperbolic umbilic meta- 
morphoses, i.e, the two symmetrical occurrences on the 
wavefront of the exchange of a swallowtail point from 
one cuspidal ridge to another. Both metamorphosis are 
part of Arnold's list of metamorphoses of fronts in space 
— see, for instance, the first and fourth perestroikas in 
Figure 258 on page 489 of @. 



2. Caustic Sheet 

If we think of the family of wavefronts for all times 
T as spatial slices of the observer's lightcone, then the 
collection of all the caustic singularities (i.e., cuspidal 
ridges, swallowtail points) on the instantaneous wave- 
fronts forms a two-surface in spacetime. The projection 
of this two-surface into the comoving space E is referred 
to as a comoving caustic surface. 

As a surface in the three-space E, the caustic surface is 
traced in time by the caustics of the traveling wavefront 
as the front moves away from the observer. The meta- 
morphoses of the wavefront 's caustics thus build up a 
picture of the caustic surface's singularities. The generic 
singularities of a caustic surface in three-space are of five 
different types: folds, cuspidal ridges, swallowtail points, 
elliptic umbilic points and hyperbolic umbilic points pl[ . 
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We showed that the caustic surface on an observer's 
lightcone is given by equation (|59|). In the particular case 
of the elliptical potential, the resulting caustic surface is 
shown in Figure 0. A significant feature of our caustic 
surface in three-space is that it consists of two separate 
intersecting sheets, one for each non-trivial root (57]) of 
the Jacobian determinant in (|5^). The two bottom pan- 
els in the figure show the two component sheets. One 
sheet starts closer to the observer as a horizontal "beak" 
and then develops a diamond cross section. The second 
sheet starts inside the first one as a vertical "beak" and 
eventually opens up into a smooth surface with an oval 
cross section. The transitions of both sheets take place at 
the points where they intersect. The points of intersec- 
tion where the vertical beak disappears from the second 
component sheet and the vertical cusp ridge appears on 
the first component sheet are hyperbolic umbilic points 
— there are two such points. 

We can also at this point give the exact value of the 
time of the first singular wavefront. This is the smaller 
of T(s+(0),z+(0)) and T(s_(0), z_(0)) since rays pass- 
ing through the origin have the longest delay (as is 
evidenced from the spike). For the previous subsub- 
scction, we have T(s + (0), z+(0)) = T ; + 9699.27 and 
T(s_(0),z_(0)) = Ti + 10148.84, i.e., the first singular 
wavefront is the surface at T = TJ + 9699.27, as antici- 
pated. 

By slicing the caustic sheet with constant-^ surfaces we 
obtain a series of planar caustics representing the meta- 
morphosis of the hyperbolic umbilic. These are shown in 
Figure ^ and coincide with the well-known caustic curves 
in the gravitational lensing literature for an elliptical po- 
tential — e.g., p. 386 of [||. Beautiful photographs of 
caustic metamorphoses are shown in Jl4| — see p. 68 for 
the hyperbolic umbilic metamorphosis. 



B. Singular Elliptical Potential Lens 

We shall illustrate the wavefront singularities and caus- 
tic surface due to a singular elliptical potential, namely, 
the case of vanishing core radius (r c = 0). In this case, 
it is simpler to use polar coordinates (r, i?) on the lens 
plane. We have for the potential, 



= A a ry/i + ecos(2tf) 
and for the associated bending angle, 
A 



• Sill a 
X 

ep 



x/l + ecos(2?9) 



(1+e) cosi?,(l-e)sintf), 



which has no dependence on r. The potential vanishes at 
the origin r = 0. This case is singular because the sur- 
face mass density associated with the potential blows up 
at the origin of the lens plane (i.e., V%^>fp 9 {r) — > oo as 
r — > 0). The center of the lens (i.e., the origin) thus acts 



as an obstruction to the light rays — see p2| ) and p. 544 
of 0] for a detailed treatment of obstruction points. Soon 
after going through the lens plane, the wavefronts evolve 
markedly differently from the nonsingular case. For our 
illustrations, we have used the following values of the 
lens parameters: Ao — 0.3, e = 0.02, zi = 4000. As in the 
regular case, these values are chosen for pedagogical rea- 
sons, although taking care to respect the approximations 
of small angles. 



1. Wavefront Singularities 

Due to the obstruction at the center of the lens (i.e., the 
origin), all the wavefronts are singular and there are only 
two regimes. The first regime consists entirely of a single 
wavefront (the earliest one to make it past the lens plane) 
with an interior spike pointing towards the observer, as 
did the wavefronts of the nonsingular elliptical potential. 
However, in this singular case the spike is conical, and 
touches the lens plane at the origin. The point at the tip 
of the spike is removed. This wavefront, T = Tj = 4000, 
is shown in Figure ^. The figure shows both a field view 
of the wavefront, and a side view of the tip of the spike 
magnified to make the conical character apparent. 

At all times afterwards, the wavefronts resemble a gob- 
let — the second regime. However, these goblets are 
markedly different from those that resulted at late times 
in the nonsingular elliptical potential case. There is a 
single diamond-shaped cuspidal ridge, which lies at the 
throat of the goblet. The rim of the goblet is removed 
and acts as a boundary of the wavefront. These goblets 
lack the interior concave surface facing the observer that 
characterizes the goblets in the nonsingular case. The 
throat of the goblet has three sheets, each typically giv- 
ing rise to an image of the source. One such wavefront 
in this regime (T = Tj+ 1000) is illustrated in Figure [l(], 
where a field view evidences the resemblance to a goblet. 
The bottom panel in the figure shows a magnified view of 
the throat of the "goblet" where the diamond structure 
is evident. 

Two questions relating to this imperfect "goblet" wave- 
front arise. In the first place, do the wavefronts ever de- 
tach from the lens plane, or, on the contrary, do they 
reach out to the obstruction point? As one can appre- 
ciate in Figure [H], where a given portion of the initial 
wavefront at T = Ti is followed up and plotted again at 
T = Ti + 1000, the wavefront does not remain attached 
to the lens plane. 

The second question is whether the rim of the gob- 
let is planar, namely, whether it lies on a constant-z 
plane. The answer is no. The rim of the goblet is a 
space curve, which is plotted in Figure |l^. It must be 
kept in mind that the rim is the boundary of the wave- 
front. This rim is removed because it is the image of the 
point r = on the lens plane under the wavefront lensing 
map, i.e., w^(0, i?) = (s(0, i?, T), z(0, tf, T)), which is a 
space curve for fixed T. Explicitly, equations (|32|) and 
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( |33| ) yield that 
«(O,0,T) = 

z(O,0,T) = 



(66) 



T + zi |&^(t?)| a /2) a£"(0) 
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(65) 



As is seen from these equations, the ultimate reason why 
the wavefronts have a whole elliptical curve's worth of 
boundary points instead of only one point is that the 
bending angle vector field 6l s ™ 9 does not depend on r. 

In rectangular coordinates, the map ol s ™ 9 is ill-defined 
at the origin r = 0. It has a singularity of the type 0/0 
since 



a, 



'(r) 



1- 



((l- e ) ri ,(l + e )r 2 ) 



One can approach the singular point from different direc- 
tions in the lens plane and hope that the limit would be 
the same in all directions. However, the limiting value of 
the a™ 9 (r) as r — > is direction-dependent, as signaled 
by the direction-dependent bending angle 1) in (pa). In 
this atypically singular case, obstructing one light ray 
(the one that passes through the center of the lens) re- 
moves an elliptical curve's worth of points on the wave- 
front. If a source is outside the rim, then one lensed 
image is seen of the source, while two images are seen if 
the source is inside the rim. This phenomenon seems to 
have been first noted on p. 188 of D for the case of a 
singular isothermal sphere (i.e., r c = e — 0), where the 



rim is a circle, as discussed in Subsection VI C 



2. Caustic Sheet 

The caustic sheet in this case is obtained by the same 
method as in the case of a nonsingular elliptical poten- 
tial. The use of polar coordinates (r, $) simplifies the 
calculation. We have 

det^fi = cos(2tf) detfj^ 
or o(r, v) 

roA-\ ( dx dy 9xdy\ 
= COS(20) U^~^J' 

where x and y are the cartesian components of the lensing 
map s z . In this case, the components reduce to 

( z (z - zi)A (l + e)\ 
— r , I cost*, 



zi v/1 + ecos(2i?) 

' z_ r _ (z - zi)A (l - e) 
zi ^1 + ecos(2?9) 



sin 1?, 



and the vanishing of the Jacobian determinant 
equation (|56"|) — becomes 



see 
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The unphysical root z — 0, which defines the plane of 
the observer, is responsible for the absence of a second 
component caustic sheet, the most notable aspect of the 
caustic sheet in the singular case. The (one component) 
caustic sheet is obtained from the non-trivial root, i.e., 



= z c (r,0) = - 



Zl A (l - e 2 ) 



(1 + ecos(2tf))2 r/zi - A (l - e 2 ) 



(67) 

A plot of the caustic sheet is shown in Figure |13|. The 
surface has four cuspidal ridges with a conical profile, 
which can clearly be appreciated in the bottom panel of 
the figure. The tip of the caustic surface lies at the cen- 
ter of the lens (origin) and is removed. Slicing the caus- 
tic sheet with constant— z planes we obtain the z-planar 
caustic curves. The planar caustics are astroid-shaped 
curves with four cusps, as can be seen in Figure [l4| A 
source inside the caustic curve has four lensed images, 
while a source outside has two. 



3. Caustic Surfing 

The singular elliptical potential provides an excellent 
opportunity to illustrate the caustic surfing scheme. We 
start by pointing out that the caustic sheet, which is 
parametrically given by ( |66| ) and (|67]), can equivalently 
be given by the following expression: 

5 (x) = ((1 + e)x 2 )* + ((l - e)y 2 Y -(2eA (z - z ; )) f = 

This expression assumes that the apex of the lightcone 
where the caustic lies is at the origin of coordinates, and 
the lens plane is at a distance zi from it. 

Suppose now that a source is moving on the light 
source plane with an approximately constant velocity 
V = (Vi, V2,0), carrying its own lightcone with it. Then 
the caustic on the source's lightcone moves as well, and at 
any given instant of conformal time t, assuming the mo- 
tion is sufficiently slow, the caustic surface can be given 
as S(x, t) = <So(x — Vt). Assume that a space-borne tele- 
scope observes at time t = an image of the source at 
peak magnification. ^From the location of the image, by 
the (comoving) lens map we can determine the location s 
of the source on the source plane. But for our purposes we 
need to switch the point of view and think of the lightcone 
of the source instead: The source lies at the origin and the 
telescope at a location (x ,y ,z ) = (-sl u / t sU l t i+ l sl ), 
where la and i s i are the distance between the telescope 
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and the lens plane, and between the lens plane and the 
source plane, resp. Because (xq, y$, z Q ) lies on the co- 
moving future caustic sheet of the source at t = 0, we 
have 

So(x ) = = ((1 + e)xlY + ((1 - e)yiy - (2eA £ tl )% , 

where we have made the substitution zi — £ t i- Now, we 
calculate the speed v mi „ that the telescope needs to stay 
on the caustic sheet during the time increment St. ^From 
our discussion in Section [y], specifically Eq. ( |63| ) , we have 

Vmin 2 2 2 x 

[(1 + e)/x ]* + [(1 - e)/y }* + [^Al/£ tl Y 

In this expression, all the symbols are known in principle. 
The values of xq and yo can be found via the lens mapping 
at the observation event, as explained above. Clearly, 
this may be a technically complex procedure, considering 
that it is only a first step in a iterative scheme to place 
the telescope on the caustic sheet. Our main purpose is 
to show that the calculation is feasible in principle. 



C. Singular Isothermal Sphere Lens 

For a singular isothermal sphere lens, we have r c = 
e = 0. The potential and bending angle vector are given 
respectively as follows: 

#£(r,0) = A r, d™(r,0) = A (cost?, sintf). 

The wavefronts are similar to those of the singular el- 
liptical potential, except that the throat of the goblet 
is a point and the singular rim of the goblet is a circle 
- see Figure [l5| The caustic surface collapses to a line 
coinciding with the portion of the optical axis beyond 
the lens plane (i.e., z > z{). Note the origin on the lens 
plane, which is a singularity of the lens, is not a part of 
the caustic line. We have that each z-planar caustic is a 
point. 



VII. CONCLUDING REMARKS AND 
OUTLOOK 

We have shown how to construct a lensing map that 
takes points on the lens plane into points on a wavefront 
surface, as opposed to a source plane. This represents 
a sort of "instantaneous" lens map, carrying a sense of 
constant time. By contrast, the standard lens map car- 
ries no sense of time at all. Additionally, the wavefront 
lensing map has an asymmetric Jacobian matrix. No- 
tice that, barring multiple sheets, the wavefront surface 



lies very close to a plane in the weak-field case, as illus- 
trated in most of our figures. In the figures the optical 
axis is magnified several times in order to appreciate the 
distance between the different sheets in the folded wave- 
front. Our intial motivation for explicitly constructing 
a wavefront- lens mapping was inspired by several works 
(e.g., H [D| |l5|, [l^, |23|). In this respect, it appears that 
an extension of our constructions in Sees. Ill and IV 
beyond the weak field domain is feasible |24], p5[ . 

We took advantage of the conformally flat nature of 
the flat FL universe in order to make use of a comov- 
ing lens equation, which in addition to being comoving 
is also a conformal lens equation. A conformally related 
flat spacetime exists, of course, in all three types of FL 
universes, so in principle, our wavefront map could be 
adapted for interpretation in the open or closed universes. 
However, in such universes the conformal factor relating 
the FL universes to the corresponding flat spacetime de- 
pends on the space point, as well as the time, and the 
translation of our conformal present proper length func- 
tion to a cosmological present proper length is not at all 
as direct as we found it in this work. Some of the sub- 
tleties that would be involved in the translation in the 
open and closed cases are treated in detail by Frittelli, 
Kling and Newman p6t , where the conformally flat lens 
map and time delay are transformed into the cosmologi- 
cal lensing map and time delay. Nonetheless, we do not 
need to rely on the conformally related flat space to cal- 
culate our present proper length function. 

The individual wavefronts for the potentials illustrated 
proved useful in visualizing the relationship between the 
location of the source in reference to the caustic, and the 
number of images observed. Additionally, the wavefronts 
of the singular potentials helped explain the anomalous 
counting of images observed by Petters, Levine, and 
Wambsganss ||, p. 188. The latter showed that in the 
case of the singular potentials there is a simple closed 
curve, that is not a caustic, but that separates regions in 
the source plane where the number of images differs by 
one, rather than two. We have here shown that such a 
curve is the boundary of the wavefront, it is not a caustic, 
and the number of images differs by one less than in the 
regular case because one whole sheet of the wavefront is 
missing due to an obstruction in the lens plane. 

Lastly, we have taken a step towards a preliminary 
scheme for caustic surfing in a meaningful and consistent 
manner. The form of the subsequent iterations remains 
an open problem, as does the implementation of the pro- 
cedure, particularly the measurement of the transverse 
velocity of the the moving source that one intends to fol- 
low. Clearly the caustic surfing proposal depends only 
on the caustic surface, and not on the method used to 
obtain it. But the point is that optimal caustic surfing 
(with the least effort) is achieved only by allowing the 
telescope to surf the caustic sheet, rather than the planar 
caustics at fixed distance from the source. One could also 
image scenarios where a telescope may not ride a caustic 
sheet, but move so as to stay inside certain chambers of 
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the caustic surface (compare with Gaudi and Gould p7J), 
possibly near higher order singularities. Future studies 
of the aforementioned issues are clearly warranted. 
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FIG. 1: A portion of the wavefront of an elliptical potential 
in the regular regime, i.e., early times past the lens plane (in 
this case, T — T t + 8900). The observer is to the right at 
(0,0,0). 




FIG. 2: The single lips regime (T = T ; + 9800). The bottom 
panel shows a front view of the lips singularity. 
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FIG. 3: Typical wavefront in the criss-cross lip regime (T = 
Ti + 10200), where two nested perpendicular cuspidal ridges 
are present. The middle panel is a vertical slice of the top 
panel, while the bottom panel is a horizontal slice. Both slices 
are through the optical axis. 



FIG. 4: The critical wavefront (T = Ti + 10420). 
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FIG. 5: Typical wavcfront early in the late regime (T = 
Ti + 10900). The bottom panel shows a front view. 



FIG. 6: Typical "goblet" wavefront in the very late regime 
(in this case, T — Ti + 55000). The middle-top panels shows a 
magnification of the throat. The middle-bottom and bottom 
panels show vertical and horizontal slices, respectively. 
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FIG. 7: Caustic sheet of the elliptical potential. The two 
component sheets are represented separately in the middle 
and bottom panels. 



FIG. 8: The metamorphosis of the planar caustics arising 
from z-plane slices of the caustic surface. 




FIG. 9: The unique wavefront in the early regime (T = T;) 
in the singular case. The bottom panel shows a side view of 
the tip of the spike. 



FIG. 10: Typical wavefront in the late regime (T = T; + 1000) 
in the singular case. The bottom panel shows a magnified 
view of the "throat" of the goblet. 
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FIG. 11: Progression of a portion of the wavefront in the 
singular case. The observer is on the right. 
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FIG. 12: The boundary of the wavefront in the critical case. 
Middle and bottom panels show side and front views. 




FIG. 14: Typical planar caustic in the singular case. 



FIG. 15: Wavefront due to a singular isothermal sphere. The 
top panel shows a field view, and the bottom panel shows a 
slice or profile of the wavefront. Notice that the wavefront is 
a surface of revolution around the optical axis. 



